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The theory of a jamming transition is proposed for the homoge- 
neous car-fohowing model within the framework of Lorenz scheme. 
We represent a jamming transition as a result of the spontaneous 
deviations of headway and velocity that is caused by the accelera- 
tion/braking rate to be higher than the critical value. The station- 
ary values of headway and velocity deviations, and time of accel- 
eration/braking are derived as functions of control parameter (time 
needed for car to take the characteristic velocity). 
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1 Introduction 

In recent years considerable study has been given to the traffic problems 
[0]. It is shown, in particular, that the jamming transition is similar to 
the conventional gas-liquid phase transition, where the freely moving traffic 
and the jammed traffic correspond to the gas and liquid phases respectively. 
The transition between them is caused by growth of car density above a 
critical value. The congested traffic flow with unstable uniform part leads 
to the formation of traffic jams where the freely moving traffic and jammed 
traffic coexists. Within the framework of Ref. p|, the jamming transition 



is represented as a first-order phase transition, whose behavior is defined by 
headway (car density) that acts as the volume (density) and by the inverted 
delay time (sensitivity parameter) that reduces to temperature. 

Our approach is to take into consideration the complete set of freedom 
degrees as equivalent variables. We obtain the self-consistent analytical de- 
scription of the jamming transition as result of the self-organization caused 
by the positive feedback of the headway deviation and acceleration/braking 
time - on the one hand, as well as the negative feedback of the deviations of 
headway and velocity - on the other one. 

The paper is organized as follows. In Section 2 the self-consistent Lorenz 
system of the governed equations for the headway and velocity deviations 
as well as for the acceleration/braking time is obtained. The jamming tran- 
sition is shown to be supercritical in character (has the second order) if a 
relaxation time for the first of pointed out quantities does not depend on its 
value; it transforms to subcritical regime with this dependence appearance. 
Section 3 deals with the determination of steady-state values for the headway 
deviation and the acceleration/braking time within adiabatic approximation. 
Out of the latter limit, the time dependencies for the headway and velocity 
deviations are studied on the basis of the phase-portrait method. Section 4 
contains a short discussion of used assumptions. 

2 Basic equations 

Within framework of the simplest car-following model the acceleration V 
of a given vehicle as a function of its distance Ax to the front vehicle is 
defined by equality V—[voptiA.x) — V]/t, where VoptiAx) — Ax/to is the 
optimal velocity function (to being a characteristic time interval), h = Vto 
is the optimal headway and r is the time of acceleration/braking needed for 
car to take the optimal velocity. It is convenient to introduce deviations 
Tj = Ax — h and v = Ax — h/to + V of headway Ax and its velocity Ax 
from corresponding optimal values h and h/to — V. Then, the flow of cars 
can be described in terms of the pointed out quantities r], v, and r. The key 
point of our approach is that the above degrees of freedom are assumed to 
be of dissipative type, so that, when they arc not coupled, their relaxation to 
the steady state is governed by the Debye-type equations with corresponding 
relaxation times tr,, t^, tr- Within the simplest approach, equations for the 
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time dependencies ri{t), v{t), and rit) are supposed to coincide formally with 
the Lorenz system that describes the self-organization process 0. 
The first of the stated equations has the form 

V=-V/tr,+V, (1) 

where the dot stands for a derivative with respect to time t. The first term in 
the right-hand side describes the Debye relaxation during time t^, the second 
one is the usual addition. In a stationary state, when r] = 0, solution of 
Eq. (|I]) defines conventional linear relationship rj = t^v, so that the headway 
deviation is proportional to the velocity deviation. 

The equation for the rate of quantity v variation is supposed to have the 
nonlinear form 

V = -v/U + QvTlT: (2) 

where are positive constants. As in Eq. ([l|), the first term in the right- 

hand side of Eq. (|^) describes the relaxation process of velocity deviation v 
to the stationary value f = determined by a time t^. The second term 
describes the positive feedback of the headway deviation rj and the time r of 
acceleration/braking on the velocity deviation v that results in the increase 
of value V and, thus, causes the self-organization process. 
The kinetic equation for the acceleration/braking time r 

f = {To-T)/tr- grVV (3) 

differs from Eqs. (|I]), (0) as follows: the relaxation of quantity r occurs 
not to the zero but to the finite value tq, representing the stationary time 
needed for car to take the characteristic velocity (in other words, Tq is the car 
characteristic); tj. is a corresponding relaxation time. In Eq. the negative 
feedback of the quantities rj and f on r is introduced to imply the decrease 
of acceleration/braking time r with the growth of the headway and velocity 
deviations {gr > is a corresponding constant). 

The equations (|I]), (||), (0) constitute the basis for self-consistent descrip- 
tion of the car- following model driven by the control parameter tq. The 
distinguishing feature of these equations is that nonlinear terms that enter 
Eqs. (H), dl) are of opposite signs, while Eq. (|l]) is linear. Physically, the 
latter means just that the velocity deviation is the derivative of headway 
deviation with respect to time. The negative sign of the last term in Eq. (|^) 



3 



can be regarded as a manifestation of Le Chatelier principle, i.e. since an de- 
crease in the acceleration/braking time promote to the formation of a stable 
car flow, the headway and velocity deviations rj and v tend to impede the 
growth of the acceleration/braking time and, as a consequence, the jamming. 
The positive feedback of r) and r on u in Eq. plays an important part 
in the problem. As we will see later, it is precisely the reason behind the 
self-organization that brings about the traffic jam. 

To explain the relaxation transition to the stable jamming state, we will 
show further that it is quite enough to use the adiabatic approximation: 
tv = 0,tr = 0. Therefore we could proceed not from Eqs. (^, (|) but from 
much simple expressions 

V = a^T]T, a„ = t^g^; r = tq - a^r^t;, = Ugr, (4) 

which are related to the stationary case i) = 0, f = in Eqs. (H), (|^) 
respectively. The equalities (§) have absolutely clear physical meaning: the 
increase of the headway deviation t] or acceleration/braking time r leads to 
growth of the velocity deviation v, whereas the increase of the headway r) 
and velocity v deviations should cause the decrease of acceleration/braking 
time r in comparison with characteristic time tq if the car flow is not broken. 
After introducing the suitable scales for quantities rj, v, t: 

Eqs. (H), ® can be rewritten in the simplest form of the well-known 
Lorenz system: 

T] = —r] + v, (6) 
e V = —v + T]T, (7) 
St = (ro - r) - r]v, (8) 

where the relaxation times ratios e = tv/tj^, 6 = t^-lt^ are introduced and 
the dot now stands for the derivative with respect to the dimensionless time 
t/Trj. In general, the system (|) - (H) can not be solved analytically, but 
in the simplest case e -C land 5 -C 1, the left-hand sides of Eqs. (|^, (|^) 
can be neglected. Then, the adiabatic approximation can be used to express 
the velocity deviation v and the acceleration/braking time r in the form of 
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the equalities (^). As a result, the dependencies of r and v on the headway 
deviation t] are given by 



Tor? 



(9) 



r 



■ 



V = 



1 + 772- 



1+V 



Note that, under t] is in the physically meaningful range between and 
1, the acceleration/braking time is a monotonically decreasing function of ry, 
whereas the velocity deviation v increases with 77 (at > 1 we have df /d?7 < 
that has no physical meaning). 

Substituting second equality (H) into Eq. (H) yields the Landau-Khalatnikov 
relation: 



with the effective potential given by 



For To < 1, the //-dependence of $ is monotonically increasing and the only 
stationary value of t] equals zero, rje = 0, so that there is no headway devia- 
tions in this case. If the parameter tq exceeds the critical value, = 1, the 
effective potential assumes the minimum with non-zero steady state headway 
deviation rje = ^/tq — 1 and the acceleration/braking time Xe = 1. 

The above scenario represents supercritical regime of the traffic jam for- 
mation and corresponds to the second-order phase transition. The latter 
can be easily seen from the expansion of the effective potential in power 
series of 77^ <^ 1: 



So the critical exponents are identical to those obtained within the framework 
of the mean-field theory . 

The drawback of the outlined approach is that it fails to account for the 
subcritical regime of the self-organization that is the reason for the appear- 
ance of traffic jam and analogous to the first-order phase transition, rather 
than the second-order one. So one has to modify the above theory by taking 
the assumption that the effective relaxation time tj^{ri) increases with head- 
way deviation rj from initial value tj^/{l +m) fixed by a parameter m > to 




(11) 




(12) 
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the final one [^]. The simplest two-parameter approximation is as follows: 

+ (13) 



tviv) 1 + iv/voY 

where < r^o < 1. The expression for the effective potential (|ll]) then 
changes by adding the term 

and the stationary values of t] are 

VT = mo[lT[l + r]lvmiro-r,)] | , (15) 
2^00 = (^0 - 1) - Tcvl, Tc = 1 + m. 



The upper sign in the right-hand side of Eq. (|T^) is for the value at the 
unstable state r]^ where the effective potential $ + A$ has the maximum, 
the lower one corresponds to the stable state r/g- The corresponding value of 
the stationary acceleration/braking time 



= : ^2 ' '"^ (16) 



l + ^oo + \/(l±Z^'-(l-%')^o 
smoothly increases from the value 



r.n = l + ^o^^ (17) 

at the parameter tq = Tco with 

Teo = (l - vl) rl (18) 
to the marginal value Tc = 1 + m at Tq = r^. 

3 Results 

The TQ-dependencies of r^e, Vj^-, and Tg are depicted in Fig. 1. As is shown in 
Fig. la, under the adiabatic condition trit^ -C t^, is met and the parameter 
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To slowly increases being below Tc, no traffic jam can form. At the point 
To = Tc the stationary headway deviation rje jumps upward to the value 
V2r]oo and its further smooth increase is determined by Eq. (|15|). If the 
parameter tq then goes downward, the headway deviation rje continuously 
decreases up to the point, where tq = Tco and rje = r/oo- At this point the 
headway deviation jump-like goes down to zero. Referring to Fig. lb, the 
stationary acceleration/braking time Te shows a linear increase from to Tc 
with the parameter tq being in the same interval. Then, after the jump 
down to the value (1 — ?7o)~^ at tq = Tc, the stationary time Te smoothly 
decays to 1 at tq ^ Tc. Under the parameter tq then decreases from above 
Tc down to Tco the acceleration/braking time t^ grows. When the point 
(|TB|) is reached, the traffic becomes freely moving, so that the stationary 
acceleration/braking time undergoes the jump from the value (|l^ up to the 
one defined by Eq. (|r^). For tq < Tcq again the parameter Tg does not differ 
from Tq. Note that this subcritical regime is realized provided the parameter 
m, that enters the dispersion law ([13|) , is greater than value 

2 

rrimin = (19) 

Clearly, according to the picture described, the jamming generation is 
characterized by the well pronounced hysteresis: the cars initially being at 
motion with optimal headway between them, begin to deviate only if the 
acceleration/braking time tq of cars exceeds its limiting value Tc = 1 + m, 
whereas the acceleration/braking time Tcq needed for uniform car flow is less 
than Tc (see Eqs. ([17|), ([ISD ). This is the case in the limit tr/^?? — > and 
the hysteresis loop shrinks with the growth of the adiabaticity parameter 
6 = tr/tr,. In addition to the smallness of 5, the adiabatic approximation 
implies the ratio t^/t,, = e is also small. In contrast to the former, the latter 
does not seem to be realistic for the system under consideration, where, in 
general, t^. So it is of interest to study to what extent the finite value 
of e could change the results. 

Owing to the condition 5^1, Eq. (|) is still algebraic and r can be 
expressed in terms of rj and v. As a result, we derive the system of two 
nonlinear differential equations that can be studied by the phase portrait 
method |p. The phase portraits for various values of e are displayed in Fig. 2, 
where the center O represents the stationary state and the saddle point 5* 
is related to the maximum of the effective potential. As is seen from the 
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figure, independently of e, tliere is tlie universal section - the "mainstream", 
that attracts all phase trajectories and its structure is appeared to be almost 
insensitive to changes in e. Analysis of time dependencies v(t) and rj{t) 
reveals that the headway and velocity deviations slow down appreciably on 
this section in comparison to the rest parts of trajectories that are almost 
rectilinear (it is not difficult to see that this effect is caused by the smallness 
of 6). Since the most of time the system is in vicinity of the "mainstream", 
we arrive at the conclusion that finite values of e do not affect qualitatively 
the above results obtained in the adiabatic approximation. 

4 Discussion 

According to the above consideration, the simplest picture of the dissipa- 
tive dynamic of traffic flow in a homogeneous car-following model can be 
represented within the framework of Lorenz model, where the headway rj 
and velocity v deviations play a role of an order parameter and its con- 
jugate field, respectively, and the acceleration/braking time r is a control 
parameter. The model is examined to show that a jam creates if the car 
characteristic tq is larger than the critical magnitude Tc- The above pointed 
out dissipative regime is inherent in the systems with small values of the 
relaxation time for acceleration/braking, being apparently a characteristic 
of a car-driver, and large ones t^, for the headway and velocity deviations. 
According to Ref. in the opposite case > t^, t^, the system behaves in 
auto-oscillation or stochastic manners. 

It is worth-while to note that the above synergetic scheme allows us to 
explain the collective phenomena of jamming transition in the A^-body prob- 
lem with N ^ oo. Then, the question arises: why exactly three variables 
(the headway and velocity deviations 77, f, and acceleration/braking time r) 
permit to explain the nontrivial behaviour of A^-body problem? The answer 
to this question gives the theorem by Ruelle and Takens: nontrivial collec- 
tive behaviour of many-body system (type of the strange attractor) can be 
represented only in the case, when number of variables is not less than three 
The interpretation of this fact is the simplest: the first of the freedom 
degrees can be chosen as the way along the phase trajectory, the second one 
corresponds to the negative Lyapunov exponent, ensuring an attraction to 
this trajectory, the third one acts in opposite manner to give repulsion. In 
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our case of the self-organization process, the second v and third r freedom 
degrees provide the positive and negative feedbacks in Eqs. (Q), (||). 

The last question to our approach is why the just Lorenz scheme allows us 
to describes main peculiarities of the jamming transition? The answer is that 
this is the simplest approach, permitting to catch on the self-organization 
effects, just as the Landau phenomenological theory of phase transitions de- 
scribes the great variety of thermodynamical transformations in the simplest 
way 0. Let us note in this connection that the effective potential given 
by the sum of equalities (pUD , ( |T4|) acts the part of the Landau free energy. 
But the above stated synergetic scheme has the principle difference from the 
Landau-type theory because the former takes into account a feedback of 
thermostat (the velocity deviation and the acceleration/braking time) with 
subsystem under consideration (the headway deviation), whereas the latter 
does not. 

Acknowledgment 

We are grateful to anonymous reviewer for constructive criticism that pro- 
motes to improve the above scheme exposition. 

[1] Traffic and Granular Flow, edited by D.E. Wolf, M. Schreckenberg, and 
A. Bachem (World Scientific, Singapore, 1996); D. Helbing, Verkehrs- 
dynamik (Springer, Berlin, 1997). 

[2] Takashi Nagatani, Phys. Rev. E 58, 4271 (1998). 

[3] H. Haken, Synergetics. An Introduction. Nonequilibrium Phase Transi- 
tions and Self- Organization in Physics, Chemistry, and Biology, 3 rd ed. 
(Springer, Berlin, 1983). 

[4] A. Vespignani, S. Zapperi, Phys. Rev. Lett. 78, 4793 (1997); Phys. Rev. 
E 57, 6345 (1998). 

[5] A.I. Olemskoi, A.V. Khomenko, JETP 83, 1180 (1996). 

[6] A.I. Olemskoi, Theory of Structure Transformations in Non- equilibrium 
Condensed Matter (NOVA Science, N-Y, 1999). 



9 



CAPTIONS 



Fig. 1. The TQ-dependencies of the stationary values of: (a) headway 
deviations rje, rj"^', (b) acceleration/braking time Tg. The arrows indicate the 
hysteresis loop. 

Fig. 2. Phase portraits in the rj — v plane at m = 1, r/o = 0.1, tq = 1.25tc 
for: (a) e = lO'^; (b) e = 1; (c) e = 10^. 
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FIG. 2 



